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An ideal phase-preserving linear amplifier is a deterministic device that adds to an input signal 
the minimal amount of noise consistent with the constraints imposed by quantum mechanics. A 
noiseless linear amplifier takes an input coherent state to an amplified coherent state, but only works 
part of the time. Such a device is actually better than noiseless, since the output has less noise than 
the amplified noise of the input coherent state; for this reason we refer to such devices as immaculate. 
Here we bound the working probabilities of probabilistic and approximate immaculate amplifiers 
and construct theoretical models that achieve some of these bounds. Our chief conclusions are the 
following: (i) the working probability of any phase-insensitive immaculate amplifier is very small 
in the phase-plane region where the device works with high fidelity; (ii) phase-sensitive immaculate 
amplifiers that work only on coherent states sparsely distributed on a phase-plane circle centered at 
the origin can have a reasonably high working probability, 
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I. INTRODUCTION AND MOTIVATION 

Classical non-inverting amplifiers take a macroscopic 
input signal, such as a time- varying voltage, and produce 
an output signal that is a rescaled version of the input 
signal. The ratio of the input amplitude to the output 
amplitude is called the gain g of the amplifier. Classi- 
cal amplifiers are used ubiquitously, e.g., to boost sig- 
nal strength for classical communications or to increase 
the power of signals driving loud speakers. In principle, 
a classical amplifier can be noise free in the sense that 
no noise is added to the input signal. The only truly 
fundamental limit on amplification comes from quantum 
mechanics. 

The canonical quantum amplifier is called a phase- 
preserving linear quantum amplifier. It takes an input 
bosonic signal and produces a larger output signal [IH3], 
while preserving the phase. The quantum constraints on 
the operation of such a device are ultimately a conse- 
quence of unitarity and can be thought as coming from 
the prohibition on transformations that increase the dis- 
tinguishability of nonorthogonal states [H [S] . The quan- 
tum constraint on a high-gain device can be expressed as 
the requirement that the amplifier must add noise that, 
when referred to the input, is at least as big as an extra 
unit of vacuum noise. A device that achieves the minimal 
added noise is called an ideal linear amplifier. 

To understand the purpose of quantum amplifiers, it 
is instructive to look at how they are used. An illus- 
trative case involves experiments probing quantum me- 
chanics at microwave frequencies. Experimenters wish to 
measure the small amplitude and phase shifts of a field 
that is used to probe another quantum system. It turns 
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out that quantum-limited simultaneous measurements of 
both amplitude and phase shifts introduce the same ad- 
ditional unit of vacuum noise as does an ideal linear am- 
plifier [6]. Thus, in principle, measuring at the input 
or amplifying and measuring at the output both provide 
the same signal-to- noise ratio (SNR); the practical ques- 
tion becomes whether it is easier to do quantum-limited 
measurement or to do quantum- limited amplification and 
subsequent measurement at the output. The answer at 
microwave frequencies is that amplifiers operate closer to 
quantum limits. 

Recently Ralph and Lund ^ proposed a device, which 
they call a "nondeterministic noiseless linear amplifier." 
The idea behind the Ralph-Lund device is that it might 
be possible to improve the SNR in some number of tri- 
als/experiments, while the device fails in the remaining 
runs. Specifically, what Ralph and Lund proposed is a 
device that takes an input coherent state |a) to a tar- 
get coherent state \ga) with (success) probability p^ and 
fails with probability 1—p,/. Such a device is even better 
than noiseless, because when the output noise is referred 
to the input, it is smaller than the original coherent-state 
noise by a factor of 1/g^. In particular, it is better than a 
device that amplifies the input noise to the output with- 
out the addition of any noise, a device that we call a 
perfect amplifier. Because it is better than perfect, we 
call Ralph and Lund's proposal an im,maculate amplifier. 
The purpose of this paper is to analyze in detail and to 
bound the performance of immaculate linear amplifiers. 

In Sec.|lT]we review recent work on deterministic linear 
amplifiers [2, which allows us to consider on the same 
footing ideal linear amplifiers and (unphysical) perfect 
and immaculate amplifiers. We use this discussion to 
motivate the idea of nondeterministic, or probabilistic, 
versions of perfect and immaculate amplifiers, and we 
use a simple uncertainty-principle argument to bound the 
working probability of probabilistic perfect and immacu- 



late amplifiers. 



Section III reviews the relation between amplification 
and cloning, thus connecting the results in this paper to 
the literature on cloning of coherent states, and Sec. |IV| 
reviews proposals for and experimental implementations 
of immaculate linear amplifiers. 

Sections \V\ and \VJ\ are the heart of the paper, the place 
where we derive bounds on the operation of immaculate 
amplifiers. Immaculate amplifiers that produce the tar- 
get coherent state exactly, but are allowed to fail, are the 
subject of Sec. |V) they are closely related to unambigu- 
ous state discrimination [5J|S], in which one discriminates 
among a set of linearly independent states exactly, but 
can declare a failure to discriminate. We use results from 
unambiguous state discrimination to bound the working 
probability of an immaculate amplifier that amplifies M 
coherent states uniformly spaced around a circle of ra- 
dius \a\ centered at the origin of the phase plane. In 
the case of many coherent states on both the input and 
output circles, i.e, assuming M :s> g^|ap, the working 
probability is bounded by 
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This success probability decreases exponentially with M 
and goes to zero in the phase- insensitive limit M — >■ oo. 
We stress that this means that an immaculate amplifier 
that works exactly on an entire circle of input coherent 
states never works. 

For an immaculate amplifier that acts on all coherent 
states on M equally spaced spokes of a disk of any radius 
|a| > centered at the origin, the success probability is 
governed by the limiting circle of zero radius and thus is 
bounded by 



P^ < 



g2{M-l) 



(1.2) 



for any M > 2. This success probability goes to zero in 
the phase-inscnsitivc limit 7\/ — > oo. 

On a more optimistic note, we also show in Sec. [V] 
that if the M coherent states are more than about a vac- 
uum unit apart on the input circle, they can be immac- 
ulately amplified with a success probability exceeding a 
half. This suggests that practical applications of immac- 
ulate amplifiers are likely to be as amplifiers that are both 
phase sensitive and amplitude specific in that they only 
work well on a discrete set of states on a particular phase- 
plane circle. Such an amplitude-specific, phase-sensitive 
amplifier might prove useful, for example, in discriminat- 
ing the coherent states used in phase-shift keying [TIIIITT] . 

The results of Sec. |V] indicate that exact immaculate 
amplification and phase insensitivity don't go well to- 
gether. In Sec. |VI| we explore this incompatibility further 
by dropping exactness and investigating the performance 
of approximate, probabilistic immaculate amplifiers that 
are explicitly phase insensitive. We characterize such a 
device by its amplitude gain and by the radius yN / g 



of the disk, centered at the origin, over which it ampli- 
fies an input coherent \a) to the target output state \ga) 
with near unit fidelity. The high-fidelity outputs thus lie 
within a disk of radius ^/N . By finding the optimal such 
amplifier, we show that the best success probability in 
the high-fidelity input region is 
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(1.3) 



which decreases exponentially with N. We use our re- 
sults to investigate the performance of phase-insensitive 
immaculate amplifiers within the context of the signal- 
to-noise ratios for measurements of amplitude and phase 
shifts discussed above. 



Because the success probability 1.3) is so small, we 



suggest that a good performance measure for phase- 
insensitive immaculate amplifiers must include both the 
fidelity with the target output \ga) and the success prob- 
ability. A natural combination is the product of the two, 
which can be thought of as the overall probability to 
reach the target. We show that over the whole range 
of operation of the optimal phase-insensitive immacu- 
late amplifier, this probability-fidelity product is never 
better than that of the identity operation. This can be 
summarized by saying that in terms of the probability- 
fidelity product, phase-preserving immaculate amplifica- 
tion is never better than doing nothing, thus reenforcing 
our conclusion that any practical application of immac- 
ulate amplification lies in phase-sensitive amplification. 

A concluding Sec. |VII| wraps up by summarizing our 
key results and discussing avenues along which future 
research might and should proceed. 



II. PHYSICAL AND UNPHYSICAL LINEAR 
AMPLIFIERS 

A. Context 

The setting for our investigation is a signal carried by 
a single-mode field, 

E{t) = - (ae"'"* + a^e"'"*) = —= {xi cosujt + X2 sinujt) . 
2 v2 

(2.1) 

This primary mode, which we label by A, is to un- 
dergo phase-preserving linear amplification. The anni- 
hilation and creation operators, a and o^, are related to 
the Hermitian quadrature components, Xi and X2, by 
a — {xi +ix2)/V^, a^ = (xi —ix2)/V2, where [a, a^] = 1 
or, equivalently, [a;i,a;2] = i- 

The annihilation operator is a complex-amplitude op- 
erator for the field, measured in photon-number units; 
the expectation value of the field, {E{t)) = Re((a)e^*"*), 
oscillates with the amplitude and phase of (a). The 
variance of E characterizes the noise in the signal; for 
phase-insensitive noise, for which ((Aa)^) = (we use 



AO = O — (O) here and throughout), this variance is 
constant in time and given by 



2{{^Ef) = (|Aa|^ 
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(2.2) 



Here (|Aap) = ^(AaAa^ + Aa^Aa) is the symmetricaUy- 
ordered second moment of a. The inequahty foUows di- 
rectly from the uncertainty principle for the quadrature 
components, ((Aa;i)^)((Ax2)^) > 1/4. The lower bound 
is the half-quantum of zero-point (or vacuum) noise and 
is saturated if and only if the mode is in a coherent 
state I a). 

The objective of phase-preserving linear amplification 
is to increase the size of the input signal by a (real) am- 
plitude gain g, regardless of the input phase, while in- 
troducing as little noise as possible. The amplification of 
the input signal can be expressed as a transformation of 
the expected complex amplitude, 



g\au 



(2.3) 



A perfect linear amplifier would perform this feat while 
adding no noise; in the Heisenberg picture, the primary 
mode's annihilation operator, not just its expectation 
value, would transform from input to output as 



fiout — go-in ■ 



(2.4) 



The second-moment noise would be amplified by the 



power gain g^, i.e., (|Aao 



g2(|Aainr). The ampli- 



fier's output would be contaminated by the same noise 
as the input, blown up by a factor of g^ , but the amplifi- 
cation process would not add any noise to the amplified 
input noise. 

There are, however, no perfect phase-preserving linear 
amplifiers; the transformation (2.4 1 does not preserve the 



canonical commutation relation and thus violates unitar- 
ity. Physically, this is the statement that amplification 
of the primary mode requires it to be coupled to other 
physical systems, not least to provide the energy needed 
for amplification; these other systems, which can thought 
of as the amplifier's internal degrees of freedom, neces- 
sarily add noise to the output. This physical requirement 
is expressed in an input-output relation [TJ [5] , 



ga-u 



Lt 



(2.5) 



where the added-noise operator L is a property of the 
internal degrees of freedom. One usually assumes that 
(L') = so as to retain the expectation-value transfor- 



mation (2.3). Preserving the canonical commutation re- 



lation between input and output requires that 
[L,L^]=g^-l, 



(2.6) 



which implies an uncertainty principle for the added 
noise, 



(|ALp)>i(g2-l), 



(2.7) 



The amplifier must be prepared to receive any input 
in the primary mode, without having any idea what that 
input is going to be. This places the restriction that the 
primary mode and the internal degrees of freedom can- 
not be correlated before amplification. The total output 
noise is then the sum of the amplified input noise and the 
noise added by the internal degrees of freedom: 



(|Aaout|'>=5'(|Aai„|2) + (|AL|2)>g2 



The lower bound follows from the uncertainty princi- 
ples (2.2 1 and (2.7). An amplifier that achieves the lower 



bound in Eq. (2.7), thus adding the least amount of noise 



permitted by quantum mechanics, is called an ideal linear 
amplifier. 



B. Ideal, perfect, and immaculate linear amplifiers 

We can formulate a more general description of linear 
amplifiers by using the formalism developed in Ref. [3], 
where we showed that for any phase-preserving linear 
amplifier, its action on an input state p of the primary 
mode can be represented by an amplifier map 



Pout 



^£{p) = TvB[S{r)p(g>aS^ir)]. 



(2.9) 



In this expression, a is the input state of a (perhaps fic- 
titious) ancillary mode B, which has annihilation and 
creation operators b and b^ , and S{r) — e'"("''^° '' ) is the 
two-mode squeeze operator. The amplitude gain is given 
by 5 = coshr, and the noise properties of the amplifier 
are encoded in cr. The main result of Ref. [3] is that the 
amplifier map is physical, i.e., is completely positive, if 
and only if cr is a physical ancilla state. 

The P function of the output state can be written as 
a convolution of the P function of the input state with 
the Q distribution of a: 



Poutm 



dh 



g4-(/3*-ga*)/yp 



Pin (a). 
(2.10) 



We specialize for the remainder of this subsection to a 
coherent-state input |a), for which the input P function 
is a (5-distribution and the output P function is obtained 
by displacing and rescaling the Q distribution of a, 



PontiP) 



Qj-{fi*~ga- 
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(2.11) 



Moments of a calculated using the P function give nor- 
mally ordered moments of a and a^ . 

An ideal linear amplifier corresponds uniquely to the 
case where the input ancilla state is vacuum, i.e., a = 
|0)(0|, giving rise to an output P function 



Pout(/3) 



-\l3-ga\'/(g' 
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(2.12) 



The displacement of the Q distribution indicates that the 



input complex amplitude is amplified as in Eq. (2.3), and 



the rescaling of the Q distribution confirms that the total 



■ Aflo 



(symmetric) output noise is (jAaoutP) = (AaQ^^iAUout/ 
i = o2 _ 1 

2 y 2- 

We can embed the ideal-amplifier map in a sequence 
of maps for both physical and unphysical amplifiers by 
considering ancilla states of thermal form, 



1- 






1 ^/ 1 



ri=0 



\n) {n\ . 
(2.13) 



When /i,^ e [l,oo), a is a physical thermal state, with 
dimensionless inverse temperature f3 given by /i^ = (1 — 
g-/3^-i. ^2 _ 2 gives the vacuum state. When ^^ £ [0, 1), 
however, a has negative eigenvalues and thus is unphys- 
ical. When /i^ g (^,oo), the trace of a is well defined 
and equal to 1, but when //^ £ [0, 1/2], the series for the 
trace of a diverges; /i^ = 1/2 makes a the parity opera- 
tor. The amplifier maps corresponding to unphysical a 
are not completely positive and thus are unphysical ^. 
In the following, we sometimes use quotes to warn the 
reader that a might not be physical. 

The Q function for a, Qa{a) = e"l"l'/''V1'M^ is well 
behaved on the entire range /i^ G (0, cxd) and becomes a 
(5-function when ^^ = 0. The output P function is the 
Gaussian 



^out(/3,/i') 



■Kfj.'^ig^ 



(2.14) 



which has normally ordered output noise (Aa^Aa) = 
fj.^{g^ — 1) and, hence, symmetrically-ordered output 



(lAaoutP) = (AaLtAaout) + ^ = fJ-Hg^ - 1) + i 

(2.15) 

The output Q distribution is 
1 



t(/3,A^') 



7r[/i2(g2_l) + l] 

(2.16) 
We now focus on three amplifiers of interest, which 
correspond to three values of /x^: 

1. The ideal linear amplifier (physical), which corre- 
sponds to /z^ = 1 and which adds the minimal 
amount of (symmetrically ordered) noise permitted 
by quantum mechanics. 

2. The perfect linear amplifier (unphysical), fj? ~ 1/2, 
whose (symmetrically ordered) output noise con- 
sists only of the amplified input noise. 

3. The unphysical /i^ = amplifier, which we christen 
the immaculate linear amplifier^ because it is better 
than perfect, and which takes an input coherent 



P function 



W function 




Q distribution 

FIG. 1: (Color online) Ball-and-stick phase-space depictions 
of input and output noise for ideal (/i^ — 1), perfect (fi^ ~ i)i 
and immaculate {fi^ — 0) amplifiers defined by the amplifier 
map (2.91 with initial ancilla "state" ( |2.13 l. Color and fill 
conventions: solid (purple) fill is used for input noise; (red) 
fill with slanted lines for the output noise of an ideal amplifier; 
(blue) fill with dots for the output of a perfect amplifier; and 
solid (green) fill for the output of an immaculate amplifier. 
The primary- mode input is a coherent state |a) with \a\ — 1, 
and the gain is g = 4, giving the output state a mean that 
lies on a circle of radius g\a\ — 4. The input and output 
states are represented by noise circles centered at the mean 
complex amplitude (the stick) and having radius E/2\/2 (the 
ball), where E^ = (| Aap) is the variance of the complex am- 
plitude calculated from the appropriate quasidistribution; for 
the normal ordering of the P function, Ep = {Aa^Aa); for 
the symmetric ordering of the Wigner W function, E^ = 
^{{Aa^Aa) + (AaAa^)) = Ep -f |; for the antinormal order- 



ing of the Q distribution. 



^l 



(AaAat) = E^ -h |. The 



P-function depiction is the one suggested by the amplifier 
map (2.9 1: the dot (Ep — 0) for the input coherent state 



\a) is amplified by an immaculate amplifier to a dot for the 
output coherent state \ga); the output for a perfect amplifier 
has additional noise Ep = ^{g'^ ~ !)> a^nd the output for an 
ideal amplifier has additional noise Ep — g'^ — 1. The sym- 
metrically ordered moments of the Wigner W function give 
the traditional picture of amplifier noise: the input coherent 
state, represented by a circle corresponding to E^ — |, has 
its noise amplified by a perfect amplifier along the (grey) ra- 
dial lines to the circle with T,'^ = \g^; the output of an ideal 
amplifier has additional noise \{g'^ ~ l)i giving total noise 
^w = ff^ ~ |) a-nd the output of an immaculate amplifier 
has its noise reduced by |((?^ — 1) to the coherent-state value 
E|i/ = I . The antinormally ordered moments of the Husimi 
Q distribution give a picture suited to discussion of simulta- 
neous measurements of the quadrature components (see text) : 
the input coherent state, represented by a circle correspond- 
ing to Eq = 1, has its noise amplified by an ideal amplifier 
along the (grey) radial lines to a circle with Eq — g^\ the out- 
put of a perfect amplifier has less noise by |((7^ ~ l)i giving 
total noise Eq = \{g'^ + l)i and the output of an immaculate 
amplifier has its noise reduced by g^ — 1 to the coherent-state 
value Eo — 1. 



state I a) to an amplified output coherent state \ga). 
We let A denote the amplifier map ( 2.9 ) for the case 



of an immaculate linear amplifier, i.e., 
-4(|a)(a|) = \ga){ga\. 



(2.17) 



The operation of these three amplifiers can be un- 
derstood intuitively in terms of how the output noise 
arises from amplified input noise and added noise. The 
three canonical quasidistributions, the P function, the 
Wigner W function, and the Husimi Q distribution [T5] . 
with their different operator orderings, quantify the noise 
differently and thus provide three different perspectives 
on the relation between input and output noise. In Fig.[l] 
we illustrate the amplification transformations for ideal, 
perfect, and immaculate amplifiers. The transformations 
can be summarized in terms of ball-and-stick phase-space 
diagrams that depict the input and output noise as cir- 
cles of uncertainty centered at the input and output 
mean complex amplitudes. We give such diagrams for 
the normally ordered variances corresponding to input 
and output P functions, as in Eq. (12. 141), and also for 



the symmetrically ordered moments of input and output 
Wigner W functions and the antinormally ordered mo- 
ments of input and output Husimi Q distributions. 

The P-function perspective, with its normally or- 
dered moments, is matched to the immaculate amplifier 
map (2.17). The immaculate amplifier takes an input 



coherent state to an amplified coherent state; in the P- 
function depiction, it takes an input dot in the phase 
plane to an output dot, without adding any noise. All 
the output noise for a perfect or an ideal amplifier ap- 
pears to be added noise. 

The symmetrically ordered moments of the Wigner 
function give the traditional perspective on amplifier 
noise. A perfect amplifier amplifies input coherent-state 
noise without adding any noise. An ideal amplifier adds 
further noise (|ALp) — }^{g^ — 1), and an immaculate 
amplifier subtracts the same amount of noise. 

The antinormally ordered moments of the Q function 
give a picture matched to an ideal amplifier. The input 
noise of a coherent state is amplified by an ideal ampli- 
fier to produce the output noise without addition of any 
further noise. A perfect amplifier has less output noise 
by |((?^ — 1), and an immaculate amplifier less noise by 



C. Naive uncertainty-principle bounds on 
probabilistic /i^ amplifiers 

The antinormally ordered noise of the Q function has 
a physical interpretation that sheds light on the perfor- 
mance of linear amplifiers. Suppose one wishes to deter- 
mine the center of a coherent state by making simulta- 
neous measurements of the two quadrature components. 
The statistics of ideal simultaneous measurements are 
given by the Q distribution [B], so in i^ such measure- 
ments, one can determine the center with uncertainty 



(i5a;i)i„/v^ — {5x2)\n/ \/^ = Vv^! the uncertainties 
here, distinguished by a 5, are calculated from the Q dis- 
tribution, i.e., using antinormal ordering. Alternatively, 
one could amplify the coherent state with an ideal linear 
amplifier and determine the center of the output state 
with uncertainty (feOout/V'^ = ('5a;2)out/V^ = g/\/v; 
this allows one to determine the center of the input co- 
herent state with the same uncertainty as measurements 
at the input, i.e., (fai)out/ffv^ = {5x2)out/ gy/v = '^/\/v. 
The point of linear amplification is to make a signal much 
larger so it can be detected by less sensitive measure- 
ments. That it is possible to determine the input with 
exactly the same sensitivity by measuring either the in- 
put or the output is an alternative way of characterizing 
the performance of an ideal amplifier. 

It is interesting to apply this sort of thinking to the 
unphysical amplifiers with /i^ < 1 ; if one could construct 
such an amplifier, one could determine the center of an 
input coherent state with uncertainty 



{5xi 



5V^ 



('^a;2)out ^ v/M'(g'-i) + i 
gy/v gV^ 



(2.18) 



This violates the uncertainty-principle bound for any 
/i^ < 1 and thus provides another way of seeing why 
the amplifiers with /i^ < 1 are unphysical. 

A potential way to make such an amplifier physical 
is to make it nondeterministic, so that it only works 
with probability p^ . Then, since only p^v of the tri- 
als are effective, one can determine the center of the in- 
put coherent state with uncertainty {S xi) out /g^/p7i^ = 
{S X 2) out /g^/Pv^- Requiring that this uncertainty not 
best the uncertainty-principle bound. 
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jives us a bound on the working probability. 
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(2.19) 



(2.20) 



Another way to express the bound (2.201 is in terms 
of the root-probability-SNR product, y'p^SNR, where 
if a;i and X2 represent the amplitude and phase quadra- 
tures ((xi) = v^laj and (^2) = 0), the signal-to-noise 
ratio is defined as SNR = {xi)/6xi ~ {xi)/Sx2- The 
root-probability-SNR product is a measure of the resolv- 
ability of states. The uncertainty-principle bound (2.20) 



on success probability is equivalent to the requirement 
that amplification not increase this resolvability, i.e., 



p7SNRout <SNRi„ = V2|a| 



(2.21) 



The root-probability-SNR product provides the same in- 
formation as the uncertainty-principle bound, but with- 
out referring output quantities to the input. We consider 



the root-probability-SNR product again in Sec. VI 



It is worth noting that since the output state pout 
is Gaussian, its fidelity with \ga) is the inverse of the 
antinormally-ordered output variances: 



F(p'^) = (galPoutlga) = T^Qoutiga) = 



1 



(2.22) 
This gives a bound on the probabihty- fidelity product, 



pAp')F{^l^)< 



g^ 



(2.23) 



which is independent of ^^ and achieved by an ideal lin- 
ear amplifier. The probability-fidelity product can be 
regarded as the overall probability to reach the target 
state \ga). Such products appear again throughout our 
analysis. 

For the remainder of the paper, we focus on the im- 
maculate linea r amp lifier (/i^ — 0), for which the prob- 
ability bound (2.20) becomes p^ < l/g^- Our analysis 



shows that a nondeterministic immaculate linear ampli- 
fier only works with high fidelity on a portion of phase 
space, where it has considerably less chance of working 
than this bound. It thus does considerably worse than 
a deterministic linear amplifier in determining the cen- 
ter of an input coherent state. This suggests that such 
devices should not be thought of primarily as linear am- 
plifiers. They could be used, however, as probabilistic, 
approximate doners, a task that we consider now. 



III. AMPLIFIERS AND CLONING 

Exact, deterministic cloning is not allowed by quan- 
tum mechanics [11[51[TJ. For coherent states, the impos- 
sibility of exact, deterministic cloning corresponds to the 
impossibility of deterministic immaculate amplification. 
If one has M clones of a coherent state |a), they can 
be coherently combined in an M-port device to produce 
M—1 vacuum states and a single amplified coherent state 
\ga), with g — vM; running an amplified coherent state 
\ga) backwards through the same device splits that state 
into M clones. This equivalence between cloning and 
immaculate amplification is the basis for links between 
cloning and amplification (see, e.g., Refs. [U HI]); here 
we summarize the links and the terminology relevant to 
this paper P^[T7]. 

The cloning literature phrases the task of cloning in 
terms of transforming N replicas of the state to be cloned 
into some number M of identical clones; this is termed 
"iV to M" cloning and is often denoted A^ — >■ M. An am- 
plifier with amplitude gain g can be thought of as doing 
1 — > M = y/g cloning. Since exact, deterministic cloning 
is ruled out by the no-cloning theorem when M > N, one 
must drop either exactness, considering instead noisy or 
approximate cloning |18| , or determinism, considering in- 
stead probabilistic cloning. 

Consider first approximate, deterministic cloning. 
The standard measure of performance for approximate 



cloning is the fidelity F of the clones with the desired 
target state. If the clones all have the same fidelity with 
the target state, the cloning process is said to be sym- 
metric. If the fidelity of the clones is independent of the 
input state, the cloning is called universal. 

It is known [121 HH] that the optimal fidelity for cloning 
coherent states \a) to M clones that have Gaussian noise 
is achieved by using an ideal linear amplifier with gain 
g = ^/m, followed by an M-port device that splits the 
amplified state into M approximate clones, each of which 
has the marginal state pa- T he st ate pa has P function 



PaW) — g^PovLt{g(3) [see Eq. (|2.12[)], and the correspond- 



ing Q distribution is 

Qa(/3) 

The output fidelity. 



H/3-a|V(2-'l/s") 

^(2-l/g2) 



Fi^M = (a|/Oa|a) = T^Qa{a) = 



M 
2M - 1 ' 



(3.1) 



(3.2) 



is a function of the gain alone, independent of the am- 
plitude of the input state [121 HO) • This output fidelity 
limits to I as M — >■ oo. 

Suppose instead that one desires perfect clones and 
is thus willing to put aside determinism. This is called 
exact {F = 1), probabilistic cloning ^2T , and the ap- 
propriate measure of performance is the probability p^ 
that the cloning process works. In probabilistic cloning, 
one usually restricts to a finite set of input states and 
attempts to clone these states optimally. The restriction 
on input states is referred to as state- dependent cloning. 

In Sec. |Vj we consider exact, but probabilistic im- 
maculate amplification. Given the equivalence between 
immaculate amplification and exact cloning, this can 
equally well be thought of as exact, probabilistic, 1 — >■ 
M — y/g cloning of coherent states. We show that ex- 
act, probabilistic immaculate amplification of all coher- 
ent states — or even of all the coherent states on a circle 
centered at the origin of phase space — cannot occur with 
a nonzero probability of success. If, however, the input 
coherent states are restricted to a finite set equally spaced 
around a circle centered at the origin, exact immaculate 
amplification can occur with a success probability given 
by the probability of unambiguously discriminating the 
input coherent states [H E] • Once one has identified un- 
ambiguously the input state, one can do any state trans- 
formation, including making an amplified coherent state 
or making as many exact clones as one wants. Thus we 
have a recipe for making an exact, probabilistic immacu- 
late amplifier or an exact, probabilistic, state-dependent 
doner. 

In Sec. |VI[ we derive rigorous bounds on the success 
probability of an amplifier that amplifies coherent states 
near the origin immaculately with fidelity near unity, but 
has output fidelity that decreases to zero as the ampli- 
tude of the input coherent states increases. Since the 
output states do not have Gaussian noise, the connection 



to cloning is not precise, but for coherent states near the 
origin, these amphfiers can be thought of as doners that 
are approximate, probabihstic, and state dependent. 

There is some cloning literature that considers var- 
ious combinations of approximate, probabilistic, and 
state-dependent cloning. For example, Chefles and Bar- 
nett [22] interpolate between exact, probabilistic, state- 
dependent doners and approximate, deterministic don- 
ers, including both fidelity and success probability as per- 
formance measures, but only for two input states, a re- 
striction that makes their results too limited for our pur- 
poses. There is also work on cloning for a distribution 
of input coherent states |23], which derives the optimal 
average fidelity of a 1 — >■ 2 doner that acts on a Gaussian 
distribution with width A centered at the origin. As the 
width goes to zero, the average fidelity not surprisingly 
approaches unity. 



IV. PRIOR WORK ON PROBABILISTIC 
IMMACULATE AMPLIFICATION 

Ralph and Lund [7, conceived the notion of an immac- 
ulate linear amplifier and proposed a probabilistic imple- 
mentation (what they called a nondeterministic, noiseless 
linear amplifier) described by a quantum operation 



^ainp(p) = £^{P) + £ta,il{p), 



(4.1) 



where £^ is the quantum operation when the amplifier 
works and £faii, the quantum operation when it fails, de- 
scribes its fallible nature. 

Ralph and Lund and collaborators [24] suggested 
that the most straightforward incarnation of a proba- 
bilistic immaculate amplifier is to have 



£^{\a){a\) =p^\ga){ga\ 



(4.2) 



for all input coherent states, where p^ is the state- 
independent probability that the amplifier works. Since 
this makes £^ — p^A, i.e., a multiple of the map (2.17) 



for a deterministic immaculate amplifier, it is not com- 
pletely positive unless the success probability is zero. In- 
deed, quite generally, if £^ works as a linear amplifier 
with uniform success probability over the entire phase 
plane, complete positivity imposes the same restrictions 
on £^ as for a deterministic linear amplifier; in particu- 
lar, famp would be just as noisy as a deterministic ampli- 
fier, the only difference being that some of the time the 
amplifier wouldn't work at all. To make an immaculate 
amplifier physical, one must make it not just probabilis- 
tic, but also drop the idea that it can work immaculately 
over the entire phase plane with uniform success prob- 
ability. In making models of immaculate amplification, 
this is precisely what Ralph and Lund did [7] . 

For the remainder of this section, we review some of 
the theoretical proposals for and experimental realiza- 
tions of Eq. (4.1). Here implementation is interpreted 



high fidelity in a restricted region of phase space near the 
origin and with the success probability p^ depending on 
the distance of the input coherent state from the origin. 

Quantum-scissors proposal. Ralph and Lund originally 
proposed to implement Eq. ( |4.2| ) using a network of beam 
splitters, single-photon sources, and single-photon detec- 
tors, as illustrated in Fig. [2] An input coherent state 
I a) is split up equally at an iV-port splitter, each out- 
put \a/vN) is processed through a modified "quantum 
scissors" (MQS) |5S], and the outputs of the quantum 
scissors are recombined at a second iV-port splitter. Suc- 
cessful immaculate amplification requires heralding on 
the MQSs so that they work correctly and on vacuum 
detection in A^ — 1 outputs of the second splitter. These 
heralding requirements mean that quantum-scissors pro- 
posal is probabilistic, and its region of high-fidelity im- 
maculate amplification is restricted by the requirement 
that jap <C A^. Even within this phase-plane region, the 
fidelity with the target state \ga) is a function of the 
amplitude |a| of the input coherent state. 

In Ref. [2^, Jeffers tried to reduce the need to make 
A^ so large by constructing a quantum-scissors device 
that works at the two-photon level, i.e., that implements 
the truncation-and-amplification transformation \a') = 
Co \0)+ci |1)+C2 |2) + 0(|a'|3) ^ Co \0)+gci \l)+g^C2 |2). 
Though this is a nice idea, there is a catch: it requires 
lossy beam splitters or a beamtritter. Numerically it was 
shown that, for |ap = 0.1, a single two-photon device 
performs better than A^ = 3 single-photon MQSs with 
respect to the fidelity of the output with the target am- 
plified state \ga) and the success probability. No mention 
is made in either Jeffers's or Ralph and Lund's work of 
how close these implementations are to limits imposed 
by quantum theory. 

Quantum-scissors implementations. One-photon scis- 
sors devices have been implemented experimentally by 
Xiang et al. [H] and Ferreyrol et al. [?7l[^ . 

The experiment by Xiang et al. |24) used an attenu- 
ated spontaneous parametric down-conversion source to 



produce an input state pin = (1 — |q;| 



+i«ni)(ii 



as meaning that the amplifier works immaculately with 



where \a\^ G [10~^,10~^]. This state is an approxima- 
tion to a uniform mixture of coherent states of fixed am- 
plitude; the motivation for considering this input state 
state was to investigate the action of the amplifier on all 
states in the mixture simultaneously. As the value of |a| 
was so small, i.e., a = a' , only A^ = 1 quantum-scissors 
device is needed. The domain of gains used in the ex- 
periment was g £ [v2; 2]. For g ~ v3, the experimental 
data showed that the amplifier was linear over the range 
|a| e [10-3,2 X 10-2]. 

Ferreyrol et al. [271 I28j implemented quantum-scissors- 
type amplifiers with A^ = 1, input coherent states with 
\a\ e [5.5 X 10^2,1], and g e [0.25,2]. Their theoretical 
modeling and experimental results are in agreement with 
the modeling and results in Ref. [25 . The first data point 
is in the region of phase space where the device has linear 
gain. Very quickly, however, the gain decreases for input 
states with |a| > 5.5 x 10-^. Their data also show that 
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tion of a single-photon addition and subtraction device 
(M = 1), which had \a\ e [0.2,1], and g G [1.25,2]. For 
input \a\ > 0.5 the fidehty of the output state with \ga) 
dropped dramatically, and the appearance of the out- 
put Wigner function departed noticeably from the target 
Wigner function in a way we return to in Sec. |VI| The 
authors point out that an equivalent quantum-scissors 
device performs worse with respect to gain and fidelity, 
both of which decrease quicker with increasing ]a] in the 
scissors case. 



FIG. 2: Device that approximates an immaculate amplifier 
(figure based on Fig. 1 of [2j). An incident coherent state is 
split equally into A'' modes at an A-port splitter. The state 
of each mode is a coherent state \a'), where a' — a/^fN\ 
N is chosen large enough that a' — a/\/N <C 1, so that 
l'^') = jo) + "*' |1) + 0(|Q^'I^) is ^^11 approximated by its 
vacuum and one-photon pieces. Each of the N modes en- 
ters a modified "quantum scissors" (MQS) |25) . shown on the 
right, which is the heart of the amplifier. When the two de- 
tectors in the MQS get results 1,0 or 0,1, the MQS is said 
to work; a feedforward phase shift (FPS) by vr, controlled 
on one of the two outcomes, is applied to the device's out- 
put mode. The result of these manipulations is that, condi- 
tioned on the MQS working, it implements the transforma- 
tion \a') — >■ (1 -|- get a') |0) = |ffc^')trunc' ^■^■' truncation of 
the state to the vacuum-one-photon sector and change of the 
relative weights of the vacuum and one-photon contributions 
so that the one-photon weight is increased; the gain is deter- 
mined by the transmissivities and reflectivities of the beam- 
splitters in the MQS. The amplified and truncated states, 
IS'-'^Otrunc' ^■'^ recombined at a second A^-port splitter. Con- 
ditional on detecting vacuum in A'— 1 outputs of this splitter, 
the output mode is in the amplified state \gQ) in the limit 
that A — > oo. Successful immaculate amplification thus cor- 
responds to heralding on the desired outcome of all of the 
MQSs, as well as vacuum detection in the A' — 1 ports of the 
final A-port splitter. 



as the coherent-state amplitude increases, the probability 
of the amplifier's working increases, and the output state 
is increasingly distorted away from the target coherent 
state. These behaviors appear in our analysis of quantum 



limits on immaculate amplifiers in Sec. VI 

Photon addition and subtraction proposals. 
Fiurasek |29] and, separately, Marek and Filip [30 
attempt to approximate the transformation in Eq. ( |4.2[ ) 
by adding and then subtracting M photons from a low- 
amplitude coherent state. The transformation for M = 1 
is aat(|o) +a|l)) ^ a(|l) + n/2^|2)) ^ [0) +2a\\), 
which has a gain of 2. This will not act like a linear 
amplifier unless |a| < 1. Generalizing to M-photon 
addition and subtraction, the gain becomes g = M + \. 
The chief problem with this method is the experimental 
infeasibility of M-photon addition and subtraction for 
M more than a very few. 

Photon addition and subtraction implementations. Za- 
vatta et al. |31| reported an experimental implementa- 



Proposals for noise addition followed by photon sub- 
traction. To overcome the difficulties of adding M pho- 
tons, Marek and Filip [3D] suggested one could simply 
add phase-insensitive noise (random displacements on 
the phase plane) and then do M-photon subtraction. 
Intuitively this can be understood as follows: adding 
noise increases the phase space area of the state; the 
subsequent photon subtraction enhances the larger pho- 
ton numbers, producing an amplified final state that is, 
roughly speaking, squeezed in the amplitude direction. 
An explicit formula is given for the success rate as a 
function of the input coherent state, M, and the mean 
number of thermal photons added. 



Implementation of noise addition followed by photon 
subtraction. Usuga et al. |32j and Usuga |33) describe 
the preparation of a displaced thermal state which is 
intended to correspond to a coherent state with added 
thermal noise. The parameters used in their experiments 
are |a| = 0.431, g S [1,2], and M S [1,4]. For 5 > 2 
(M > 1), the authors found the probability of success 
decreased drastically, and the state started to deform. 



Discussion. From the theory and experiments summa- 
rized above, several conclusions can be drawn. First, all 
of the devices produce an output state with high fidelity 
to the target coherent state \ga) only over a restricted 
region of the phase plane centered on the origin. Second, 
although the theoretical proposals allow for high gains 
and high input amplitudes, current implementations are 
restricted to small gains .g ^ 2 and small input ampli- 
tudes |a| ^ 2 by technical limitations. Third, even for 
these small gains and small input amplitudes, these de- 
vices fail almost all of the time. 



Most previous work on this subject has focused on 
proposing and analyzing the performance of specific 
schemes for probabilistic immaculate amplification. We 
take a different tack: we provide a general analysis of 
the performance of any device that attempts to approx- 
imate immaculate linear amplification. We characterize 
the amplifier by its gain and the region of the phase plane 
over which it operates with high fidelity, and we derive 
fundamental quantum limits on the probability that the 
amplifier works. 



USD BOUNDS ON PROBABILISTIC 
IMMACULATE AMPLIFICATION 



Quantum state discrimination is a decision-theoretic 
task in which an agent, who has the abiUty to perform 
any measurement he wishes, is handed a single state 
drawn from a known set of states and is told to determine 
which of the states he received. Our chief interest here 
is unambiguous state discrimination (USD): the agent is 
told never to misidcntify the state, at the cost of sure 
and sudden death, but is allowed throw up his hands in 
despair and refuse to make a decision. A set of states 
can be discriminated unambiguously if and only if they 
are linearly independent [5] ; there is a nonzero probabil- 
ity for no decision unless the states are orthogonal. In 
this section we apply USD bounds to the performance of 
exact immaculate amplifiers. We use the USD formalism 
in two ways. 

The first is to provide upper bounds on the working 
probability of an immaculate amplifier. Let p(/) be the 
probability that an immaculate amplifier works exactly 
on a set of input coherent states. Suppose that P^ is the 
optimal probability for discriminating the input states 
and P^ is the corresponding optimal probability for dis- 
criminating the amplified states. The amplified states, 
being further apart on the phase plane than the input 
states, are easier to distinguish, so P^ > P^. The over- 
all probability of successfully discriminating the ampli- 
fied states is p{/)P^. Since P^ is optimal, the ampli- 
fication process cannot increase the distinguishability of 
the states, so we must have P^ > p{/)P^. The result is 
a strict upper bound, p(/) < P^/P^, on the probability 
that the immaculate amplifier works; we cannot warrant, 
however, that this upper bound can be achieved. 

The second way we use the USD formalism is to con- 
struct models of immaculate amplifiers that have an 
achievable working probability. Once one has used USD 
to identify one of the input states, one can perform any 
unitary transformation on that state. This procedure al- 
ways produces the right transformed state when it makes 
a decision; consequently we call it, somewhat cumber- 
somely, an exact, finite-state, probabilistic state transfor- 
mation. The transformation could be the displacement 
of a coherent state required to amplify it. Since the opti- 
mal USD discrimination probability P^ can be achieved 
in principle, the result is a model for an immaculate am- 
plifier that works with probability P^ on a finite set of 
input coherent states. We call such a model a, finite- state, 
probabilistic immaculate amplifier. 

We note this formulation and subsequent analysis is 
similar to the analysis performed by Dunjko and Ander- 
sson in Ref . [34] . Their results are not explicit about the 
dependence of the success probabilities on gain and input 
amplitude, whereas we are. 



A. Helstrom bound for two coherent states 

Before turning to USD bounds on immaculate ampli- 
fiers, we consider a related bound provided by the min- 
imal error probability in discriminating two nonorthog- 
onal states. Consider two coherent states, \a) and |/3). 
A measurement that minimizes the chance of incorrectly 
identifying the state is known as a Helstrom discrimina- 
tion measurement |351 [35] . The probability of successful 
identification is 






(5.1) 



where the superscript "B" reminds us that this proba- 
bility is before immaculate amplification. It is apparent 
that as the separation, |a — /3|, between the two states 
grows, the states become orthogonal, and the probabil- 
ity of successful discrimination approaches unity. In con- 
trast, when la — /3| — )■ 0, the success probability limits to 
guessing. 

Now we use the above-described procedure, modified 
to Helstrom discrimination, to bound the working prob- 
ability p{/) of an immaculate amplification device. The 
device takes |a) to \ga) and |/3) to \g/3). Amplification 
increases the distinguishability of the states so that the 
probability of successful identification of the state is 



P' 



^(/) = ^ (l + Vl - e-s'\--P\') , (5.2) 

where the superscript "A" reminds us this is after am- 
plification. The overall probability to identify the input 
state correctly after amplification is 

= i (l + p{/)Vl-e-9-\a-P\-^ . (5.3) 

Since the probability for successful discrimination cannot 
increase, we must have Phoi(^) — -^HoiC^)' which gives 
an upper bound on the amplifier's success probability. 



p(/)< 



\ 



l-e-|"-/3P 



(5.4) 



This bound, which holds for any pair of states, has its 
minimum value when the two coherent states become 
very close to each other, i.e., |a — /3| — >• 0; in this case the 
bound on the working probability becomes 



pHcl 



< 



1 

9 



(5.5) 



For constructing models of immaculate amplifiers, 
Helstrom-type discrimination has the problem that it 
sometimes misidentifies the input state. Such misidenti- 
fication inevitably leads to noise in the amplifier output, 
which cannot be part of a model of an exact immaculate 
amplifier. 
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B. USD bounds 



Two coherent states 



Unambiguous state discrimination does discriminate 
states without error, but this providence requires a sac- 
rifice, namely, the no-decision measurement result. For 
two input states, \a) and |/3), the probability of success- 
fully identifying them is [36] 



-fuSD(^) 



1 



(/3|«> 



1 



M-fi\^ 



(5.6) 



In this expression, as in the Helstrom case, it is ap- 
parent that as the separation, \a — /3|, between the 
two states grows, the probability of discrimination ap- 
proaches unity. When the states get close together, 
ja — /3| — >■ 0, the probability of successful discrimination 
goes to zero. 

After amplification we have a discrimination probabil- 
ity, 

pA(/) = l-|(g/3|5a)|2 = l _ e-s'l"-'5|% (5.7) 

and an overall probability for successfully identifying the 
input state. 



^usd(^) 



p{/)p\/). 



(5.8) 



Since amplification cannot increase the distinguishability 
of the states, we have -PugD(/) < -Pusd(^) ^^^ thus an 
upper bound on the working probability. 



P(/)< 



Ajsd('^) 



1 



-\a-P\' 



^-g2\a^fi\2 



(5.9) 



USD to more than two input states. Indeed, we work to- 
ward a phase-insensitive amplifier, which must act sym- 
metrically on all input coherent states with the same 
|a|. Thus what we do is to consider a set of M coherent 
states, \aj) — |ae"*j), all located on a circle of radius a 
with phases 



2^ 
M 



j=0,l,2,...,M-l, 



(5.11) 



distributed uniformly around the circle. To avoid clutter 
in what follows, we use, as here, a = \a\. To apply USD 
to the states \aj), they must be linearly independent. 
This property was shown in Ref. |37| , and it emerges nat- 
urally as part of the USD construction. In contrast, the 
continuum of states on the circle are complete, spanning 
the entire Hilbert space, but are not linearly independent; 
we review these facts in Appendix \K\ 

Chefles and Barnett [5] solved the USD problem 
for sets of linearly independent symmetric states (see 
also [8]). For the case of coherent states on a circle, 
the unitary operator that rotates between states is the 
phase-plane rotation by angle 27r/M, i.e., U = e*^'^° a/M _ 
Restricted to the subspace spanned by the set of input 
coherent states, U has the eigendecomposition 



U 



M-l 
r=0 



\lr) (7r| 



(5.12) 



where the (orthonormal) eigenstates are given by 



M-l 



Cr \^r) 



M ^ 



e-*""'^^/*^ \a. 



(5.13) 



i=o 



as pointed out in Ref. [21]. Being the square of the 
Helstrom bound (5.4), this is always the tighter bound. 



The minimum of the bound is found in the limit that 
the coherent states become very close to each other, i.e., 
ja — /3| "> 0, in which case the bound becomes 



1 

PUSD < ^, 



(5.10) 



The allowed working probability is a factor of 1/g smaller 
than the Helstrom bound (5.5 1. This USD bound is the 



same as the bound (2.201, which was derived by consider- 



ing how to distinguish neighboring coherent states using 
quadrature measurements; the two bounds are the same 
because both are based on discriminating neighboring co- 
herent states. 



2. M coherent states on a circle 



The USD bound (5.10) is not at all a tight bound on 



the working probability for a probabilistic immaculate 
amplifier. We can get much tighter bounds by applying 



Here c^, chosen to be real, is the magnitude of the vector 
on the right: 



A/-1 



cl = 



_l_^^-.2..,/M^^^|^^.^ 



J=0 



M ^ 



-ir(pj 



exp[a2(e*'^^- - 1)] 



(5.14) 



i=o 



It is useful to manipulate c^ into a quite different form 
and also to write it in terms of 



Mci 



jM-r 



dx^'-^ 



Me-' V -^ 



M-l 

y exp^xe 

2{kM+r) 



,i<t>j ^ 



fc=0 



{kM + ry. 



(5.15) 



That the states \^r) are orthonormal establishes that 
they and the original coherent states \aj) span an M- 
dimensional subspace and thus that the \aj) are linearly 
independent. 
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The vectors 



Af-l 






r]/M 



M '—' Cr 
r=0 



\lr) 



(5.16) 



are reciprocal (or dual) to the original coherent states 
in the sense that (a^jafc) = iJ^^. This duality property 
is what is needed to construct the USD POVM. This 
POVM has M POVM elements E^ = P(/) \a\) {a\\, 
J = 0, . . . , M — 1, for the results that identify the input 
states, where P(/) is the success probability, and a single 
failure POVM element, E'faii = I — E, where 



E 



^ii;,^P(/)^l|7.)(7. 



(5.17) 



The largest eigenvalue of E must be no larger than 1, 
which gives an optimal success probability for discrim- 
inating among A'l coherent states symmetrically placed 
on a circle of radius a ^ : 



P{/\a,M) 



mm Qr 
re{a,...,M-i} 



(5.18) 



This success probability has two important limits: 
(i) many states on the circle or, equivalently, small 
coherent-state amplitude, i.e., M ^ a^, and (ii) states 
sparse on the circle or, equivalently, large coherent-state 
amplitude, i.e., M <C a. The reason for the difference 
in powers of a in the two limits emerges as we examine 
each limit in turn. 



Notice first that the sums for Qr/M in Eq. (5.15) con- 
sist of terms drawn with period M from a Poisson dis- 
tribution that has mean a^, a distribution we denote 
throughout by Pr[n | a^ ] = e~" a^"/n!. When the first 
term in the sum for r ~ M — 1 lies beyond the maxi- 
mum of the Poisson distribution, as it does in the case 
of many states on the circle, it takes only a moment's 
contemplation to realize that the terms in the sum for 
qm-i are term by term smaller than the corresponding 
terms in the sums for other values of r, provided that the 
first term in qm-i is smaller than the first term in go. 
I.e., a 



2(Af-l) 



/{M — 1)! < 1, which is certainly true when 
M ^ a^. Thus, for many coherent states on the circle. 



the minimum in Eq. (5.18) is achieved by r ^ AI — 1 [S], 
so 



P(/|a,Af) = gAf-i=Me-"'^ 



fe=0 



Q,2(fcM+M-l) 

{kM + M -1)1 



(5.19) 



Moreover, the Chernoff bound for a Poisson random vari- 
able n with mean a^ [38 , applied to the terms in the 



sum (5.19) after the first, 



°°^ ^2(fcM+M-l) 



E 



^(fcM + M-l) 



< e" Pr[ n>2M ~l\a^] 



< 



\2M- 1 



(5.20) 



shows that, in the limit M ^ a^, we need to keep only 
the first term. A; = 0, of the sum (5.19). The result is a 



simple expression for USD success probability in the case 
of many coherent states on a circle (small coherent-state 
amplitude) : 



P{/\a,M) = 



(M-1)! 



M:^ a'^ 



(5.21) 



Now consider the case of sparse coherent states on the 
circle. For fixed M, as a — ?> oo, Chefles and Barnett [5] 
showed that all of the qr limit to 1, so 



P{/\a,M) = l 



(5.22) 



Since, for fixed M, the input states limit to being orthog- 
onal as a — >■ oo, this simply means that orthogonal states 
can be discriminated with unity probability of success. 
More useful than the limit, however, is the correction to 
the limit. 

To find this correction, we begin by noting that since 
a ^ M > 2, we can approximate the Poisson distribu- 



tion in Eq. (5.15) as a Gaussian of the same mean and 



variance and extend the sum on fc to — oo on the grounds 
that the Gaussian is negligible for these additional terms: 



E -P 2cP 



27ra 



(5.23) 



By introducing (5-functions, we can write this in the form 



M 



Qr = 



2TTa 



y^ /.oo 

Y^ / dxe-(^-"')'/2"'j(a;-fcM-r) 

._ J —OO 



k— — oo 



due-^"'^ Y 5[k 



/i:— — oo 



a s — u 

— H 

m m 



(5.24) 



where a; is a continuous version of kM -\- r, and where in 
the second expression we introduce the integration vari- 
able u = x/a — a and rescaled variables m = M /a <C 1 
and s = r/a <^ 1. Now we write a/m = [a/m] -\- H, 
where [z] denotes the nearest integer to z and, hence, 
— ^ < H < i (half-integers are rounded up), redefine the 
dummy summing variable to be fc — [a/?n], and use 



E ^c^ 



k=-c 



E' 



-i27TJV 



(5.25) 



to put Eq. (5.24) in the form 
1 "" 



/27r 



:[. 



OO 

2\^cos 

s 



y^ ^i27rj{s/rn-H) / ^^ ^~u^ /2 ^-i2-,rju/m 



H 



271 j[ - 
m 



-2TX^-lr. 



n 1 -1 I 2 

-271 J / m 



(5.26) 
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Here 6*3 denotes a Jacobi theta function. 

When ?7i <C 1, we only need to keep the j = 1 term 
in the sum to get the dominant correction to unity in q,,. 
To minimize g^ , we choose r/M — H = s/in — H as close 
to I as possible, consistent with letting r be an integer. 
Thus we choose r — [il/(H + |)], which gives 



cos 



2tt 



-N 



= -1 



irrelevant errors of size < 



2M2 



(5.27) 
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Keeping more terms in the sum and then minimizing 
could provide a better approximation, but the lowest- 
order, j ~ I correction already provides a good approx- 
imation for a reasonably dense set of coherent states so 
the following analysis is restricted to it. 

The resulting success probability in the case of sparse 
coherent states (large coherent-state amplitudes) is 



FIG. 3: (Color online) Dependence of the ratio o? jM^ — a(e) 
on the deviation t of the success probability P(/|a, M) from 
unity: numerical results are plotted as (red) circles; analytic 
approximation of Eq. (5.291 as (blue) squares. The analytic 



P(/|a,M) = l-e~ l-2e 



-2-K^B? jM'^ 



M <^a. 



(5.28) 



The key result here is that in this limit the success proba- 
bility only depends on the ratio a/M . Indeed, using this 
expression, we can turn the question around and deter- 
mine the ratio that gives a deviation e: 



approximation works quite well for e G [0,0.5], but breaks 
down progressively beyond e = 0.5. 



amplifiers and, second, to constructing a model of an 
immaculate amplifier based on USD. For the first task, 
we use the same notation as previously for before and 
after probabilities of USD; the USD bound on the success 
probability of an immaculate amplifier that works on the 
M input coherent states is 



p(/|a,M) < 



-PusdC^) 



a 
M2 



a{e) 



ln(e/2) 
27r2 



P(/|a,M) 
P^{/) ^ P{/\ga,M) 



(5.31) 



-0.05066 In e 



0.0351. 

(5.29) 



The important cases of this bound require only our ap- 
proximate results for the USD success probabilities. 

A first such case is when the input coherent states are 
sparse and, hence, so are the amplified output states. 
In this case, the numerator and the denominator in the 



For example, to achieve P{/\a,M) = 0.9 fo r any M, 

one chooses 6? ~ 0.15M2. The dependence (5.291 has 

been tested numerically over the ranges e e [0.5, 10"^] bound ( |5.31[ ) are both close to one, and^ the bound on 

and M G [2,40]; the numerics give 



a{e) 



-0.05081ne + 0.035. 



(5.30) 



in good agreement with the analytic approximation. Fig- 
ure [3] compares the numerics with the analytic approxi- 
mation; the analytic approximation works quite well for 
e < 0.5. 

Figure|4]plots the success probability for USD of coher- 
ent states on a circle, comparing the exact, numerically 
determined result with the approximations that apply for 
many coherent states and sparse coherent states. The 
two approximations work better than we have any right 
to expect: the plots and a consideration of the next term 
in the sum ( 5.19 1 sug gest that the many-coherent-states 
approximation (5.21) works well for M > 2a^; provided 



success probability is also close to one, reflecting the fact 
that one can discriminate and amplify such nearly or- 
thogonal states. 

More interesting is the case of many input coherent 
states, M ^ a^. If the gain is large enough that the 
amplified states are spar se, i.e ., M ^ ga — this requires 
that g ^ a — the bound (5.31 1 reduces to 



p(/|a,M)<P^SD(/) 



j^^g-a^-2(M-i) 
(Af-1)! 



(5.32) 



a is somewhat bigger than 1, the sparse-coherent-states 
approximation (5.281 works well for M < 4a. The two 



approximations overlap when a > 1 and M are both 
small, but because of the different powers of a in the two 
approximations, generally there is a gap between the two 
that must be filled in with numerics. 

These results in hand, we can apply them, first, to 
obtain bounds on the success probability of immaculate 



This bound, which is plotted in Fig. 4] as (red) circles 
in the left column and a (red) dashed line in the right 
column, can be regarded as the g -^ 00 bound on an 
immaculate amplifier that works on a fixed number M ^ 
a^ of input states. 

Most interesting is the case in which M is large enough 
that both the input and amplified output can be treated 
in the many-coherent-sta tes l imit, i.e., M ^ g^ct^- 
In this case, the bound (5.311 becomes p{-/\a,M) < 
g(s -i)a ig2{M-i)^ This case is the most interesting be- 
cause we can let M become arbitrarily large and thus 
approach the limit in which the amplifier acts phase- 
insensitively on the entire circle of coherent states. Since 
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FIG. 4: (Color online) Success probability P{/\a,M) = 
P^(/). Left column: as a function of M with fixed 6?\ 
(black) asterisks are the exact, numerically determined suc- 
cess probability ( |5.18| l; (red) circles give the approximate re- 
sult (5.21 1 for many coherent states (small coherent-state am- 
plitude]; (blue) squares give the approximate result (5.281 for 
sparse coherent states (large coherent-state amplitude). Right 
column: as a function of 6i^ with fixed Af ; (black) solid line 
is the exact result; (red) dashed line, many coherent states; 
(blue) dotted line, sparse coherent states. 



M - 1 > (g^ - l)a^, we have e'^^^''^^'^^ < e^"^ and thus 



p(/|a,M)< 



.(9^-i)a^ 



9 



2(M-1) 



< 



2(M-1) 



(5.33) 



This shows that the success probabihty of an exact im- 
maculate amplifier goes to zero in the phase-insensitive 
limit ill -^ oo, even when the amplifier is only required 
to work on a single circle of input coherent states. 



We can make a more precise statement for an immac- 
ulate amplifier that amplifies exactly all the coherent 
states on M spokes spaced equally in angle and of length 
a. Such an amplifier acts immaculately on M coherent 
states on all circles with radius < a. The success proba- 
bility is bounded by the a -^ limit of the bound ( 5.33 1, 



where the assumptions underlying the bound are satisfied 
for any M > 2: 



1 



(5.34) 



This is one of the two chief results of this section: an im- 
maculate amplifier that works exactly on M spokes within 
a phase-space disk centered at the origin has a working 
probability that decreases exponentially with M , with the 
base of the exponential, g^ , given by the gain, and goes 
to zero in the phase-insensitive limit M — > oo. 

It is useful to pause here to relate these results to 
the discussion at the end of Sec. |Tl] For the disk ampli- 
fier, the measurement-based performance measure ( 2.19 1, 



which uses antinormal ordering to calculate the uncer- 
tainties, is 1/pusd5^ > (^2(M-2). ^]-^jg jg greater than 
the uncertainty-principle lower bound of one, achieved 
by an ideal linear amplifier, for M > 2 and far worse 
than the bound as M gets large. (These same arguments 
hold for the bound on the root-probability-SNR prod- 
uct, which is equivalent to the uncertainty bound.) The 
related probability-fidelity product is given by pusD < 
l/ff^ I since an exact immaculate amplifier has unit 

output fidelity; this is worse than the probability-fidelity 
product 1/g^ achieved by an ideal linear amplifier for 
M > 2 and far worse as M gets large. 

As we discussed in the introductory paragraphs of this 
section, we can construct a USD-based model of an im- 
maculate amplifier in which the M input coherent states 
are first discriminated and then the identified input is 
amplified immaculately by any amount. The quantum 
operation for this model is 

M-l 

a(p) = J2 Pi^l^^M) \ga,) {af\p\af) {ga^l . (5.35) 

This map can be applied to any input state, not just 
the M coherent states used to construct it, but applied 
to one of those special input states, \aj), 9i outputs the 
amplified state \gaj) with probability 



p(/|a,M) = P^sd(^) = PU\a,M). 



(5.36) 



This success probability is plotted in Fig. |4j 

When M ^ a^, the success probability is given by 
Eq. (15:211), 



p(/|a,M) 




ea 
M - 1 



A/-1 



(5.37) 



where here we apply Stirling's approximation to the fac- 
torial to make clear that the success probability goes to 
zero in the phase-insensitive limit Af — >■ cx). 
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The case of sparse input states is where immaculate 
amphfication shines with the radiance its name evokes. 
As the plots in Fig. |4] show, the success probability for 
this case is captured by the sparse-states approxima- 
tion (5.28), which is plotted in Fig. H as (blue) squares 
in the left column and a (blue) dotted line in the right 
column. The approximation works well for success prob- 
abilities 1 — e ^ 0.5, which corresponds to M < Aa. 
To achieve a success probability 1 - - e re quires that 
a/M = \/a{e) be chosen as in Eq. (5.29). To get a 
feeling for what these results mean, notice that a suc- 
cess probability of 1 — e corresponds to a distance be- 
tween states, measured along the arc of the circle, given 
by 2T:a/M — 27r-y/a(e); for example, a success probabil- 
ity of 0.5 corresponds to ^/a ~ 0.265 and a distance of 
about 1.67. These states might seem pretty crowded, but 
the distance makes sense when compared with the one- 
standard-deviation diameter of a coherent state, which is 
1. These input states are just beginning to overlap, but 
they are far enough apart that they can be distinguished 
and amplified immaculately half the time. 

The lesson here is important: USD-based devices can 
outperform ideal linear amplifiers if they are both phase- 
sensitive and amplitude-specific, amplifying immaculately 
only a relatively sparse set of input coherent states on a 
particular input circle. This realization leads to a set of 
interesting questions that we consider briefly in the Con- 
clusion as the basis for future work. The flip side is that 
success probability goes to zero when an exact immacu- 
late device is required to work phase-insensitively on even 
a single input circle. This suggests that phase insensi- 
tivity is a key property, which does not play well with 
exact immaculate amplification. In the next section, we 
explore this further by considering probabilistic immacu- 
late amplifiers that are required to be phase-insensitive, 
but unlike USD-based amplifiers, are not exact. 



VI. BOUNDS ON PHASE-INSENSITIVE, 
APPROXIMATE, PROBABILISTIC 
IMMACULATE AMPLIFICATION 

In this section we canonize phase insensitivity as a pri- 
mary requirement for amplification. This means that 
the amplifier's operation must be invariant under phase- 
plane rotations. We relax the requirement of unit fidelity 
with the target output state, thus obtaining a model of 
an approximate immaculate amplifier. We would like the 
amplifier to work with high fidelity for input coherent 
states I a) within a disk centered at the origin, but we 
allow the fidelity with the target amplified state \ga) to 
fall off for inputs outside the disk of interest. There are 
two motivations for this assumption: first, as we noted 
in Sec. |IV[ an immaculate amplifier cannot work over the 
entire phase plane; second, as was true for the implemen- 



tations reviewed in Sec. IV such a cutoff is a property of 
practical devices. 

We characterize the high-fidelity output region as a 



disk of radius vN; the corresponding input disk thus 
has radius yN / g. After translating this description into 
the language of amplifier maps and Kraus operators, we 
characterize the amplifier in terms of the fidelity with 
the target state, F{a), and the probability that the am- 
plifier works, p{/\q), both of which arc functions of the 
input amplitude a = |a|. We maximize the fidelity at 
each a given a working probability at that a, after which 
we maximize the working probability consistent with the 
amplifier's map being trace decreasing. We thus obtain 
an optimal immaculate amplifier that is both approxi- 
mate and probabilistic. 

We note that a similar analysis has been performed by 
Fiurasek [40l |4T] in the context of cloning and arbitrary 
state transformations; we point out below similarities to 
and differences from our analysis. 

We describe the amplification process by a quantum 
operation, which we write in terms of a canonical Kraus 
decomposition in which the Kraus operators are orthogo- 
nal. We assume that these Kraus operators have the form 
PnKj, where P/v is the projector onto the subspace Sn 
spanned by the first A'^ -t- 1 number states. The amplifier 
quantum operation is thus 



AN = Y.Pr^^3®^]PN- 



(6.1) 



where the 0, technically a tensor product, can be re- 
garded as designating the slot for the input to the quan- 
tum operation. The projector Pn provides a sharp cutoff 
in the number basis, beyond which the amplifier's out- 
put has no support; notice that we can let the operators 
Kj map outside Sn without having any effect on the 



quantum operation ( |6.1[ ). Shortly we extend the Kraus 
operators in a way that allows the outputs to have sup- 
port outside Sn] this extension smooths the rough edges 
in the amplifier map (6.1), and it provides marginal im- 



provements in the output fidelity. Phase insensitivity is 
the requirement that An commutes with phase-plane ro- 
tations; this implies, as we show in Appendix IB] that 
each Kraus operator has nonzero number-basis matrix 
elements on only one diagonal strip, as in Eq. (B5). Ad- 



ditionally, the Kraus operators must satisfy the trace- 
decreasing requirement. 



Y,kIPnK,<I. 



(6.2) 



Suppose now that the input state to the amplifier is a 
coherent state \a). The probability of outcome j is 

p,{/\a)^{a\K}PNKj\a), (6.3) 

and the overall success probability is 

p{/\a) - 5]p,(/|«) = tr[AN{\a){a\)] . (6.4) 

j 

The fidelity of the output with the target output state 
\ga) is 



F(a) 



{ga\AN{\a){a\)\ga) 
WW) 



(6.5) 
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Because of the rotational symmetry, these quantities de- 
pend only on the magnitude a = \a\. 

The problem we solve is the following: fix a circle of 
coherent states with amplitude a, and find the maximum 
fidelity F{a) on this circle for a fixed success probability 
q = p{/\a). We do this first for a single Kraus operator 
and later argue that a single Kraus operator is better 
than more than one. The optimization problem is thus 
to maximize 



F(a) = 



\{ga\PNK\a)\^^ 



(6.6) 



Pi/\a) ' 
subject to the constraint 

q^p{/\a)^{a\K'fPNK\a). (6.7) 

We can, of course, rephrase this as maximizing 
|((?a|PAri^|a)p subject to the constraint on working 
probability. 

Introducing a Lagrange multiplier /i, we maximize 

\{ga\PNK\a)\^ - ^i{{a\K^ P^Kla) - q) . (6.8) 

Varying K gives 

^{a\6KUPN\9a){ga\PNK\a) - nPNK\a)] 

+ (Hermitian conjugate) , (6.9) 

so we conclude that 

PNK\a) = PN\ga) ^-^ ' ^ ' V (6.10) 

^^ 

The Lagrange multiplier is given by the probability for 
the first A^ + 1 photons in the target state \ga) , 

li = {ga\ Pn \ga) = e-»'l"l'ejv(ff'|an , (6.11) 

where we introduce a standard shorthand for the first 
A'^ + 1 terms in the expansion of the exponential, 



N 



Cn 



(-)-E 






(6.12) 



Without changing the Kraus operator PnK, we can let 
K map outside the subspace Sn in such a way that 

K\a) = \ga) (g"l^^^l") . (6.13) 



Since 



f*°|a)=e('^'-i)l"l'/2|^^^^ (gj4) 



we can simplify this by letting K = Lg" °. The result 

{ga\PNL\ga) 



L\ga) = \ga)- 



i" 



(6.15) 



says that \ga) is an eigenstate of L. Since the coherent 
states on a circle are a basis for the Hilbert space, this 
determines L to be a function of the annihilation operator 
a. The rotational symmetry further requires that L have 
number-state matrix elements on only one diagonal strip, 
implying that L — Aa*^, where fc is a nonnegative integer 
and A can be taken to be real without loss of generality. 
The possible optimal Kraus operators are 

Kt = Xal-g''' = Xg-'-iga)'' 



X^g'^J^^lnHn + kl 



ri=0 



fc = 0,1,2, 



(6.16) 



This operator has nonzero matrix elements only on the 
fcth diagonal strip above the main diagonal. It is not sur- 
prising that this class of operators emerges, because they 
do take \a) to a multiple of \ga), just as we would like 
an immaculate amplifier to do. The success probability 
and fidelity become 

p{/\a) - AV'e(»'-i)^'a2'^(<?a|PAr|.ga) , (6.17) 

F{a) = {ga\PN\ga)^e-s"^\N{9^a'')^^i. (6.18) 

We can increase the success probability without chang- 
ing the fidelity by letting A^ increase, but there is a limit 
to this increase set by the requirement that 

JV 



I>KlP^K,^X'y^g^(^^+'^^^^^^^±P^\n + k}{n + k\. 

(6.19) 



n=0 



Since the eigenvalues increase with n, the constraint is set 
by the largest eigenvalue {n = N). Choosing the largest 
possible value. 



A^ = 



m 



1 



(7V-f-fc)!g2(JV+fe) 



(6.20) 



maximizes the success probability. 

The final results of these considerations are the Kraus 
operators 



K.= 



{N + ky. g^+^ 



{N + k)l g^ 



(6.21) 



and the corresponding success probability and fidelity, 

TV' e""^a^'' 
P(^l«) = f.r ,'.M „2/v eA^(g« ) , (6-22) 

(6.23) 



{N + k)\ 527V 
F(a)=e-3'"'ejv(ff'«'). 



Equation ( |6.23 l was derived by Fiurasek [40] (our am- 
plitude gain g is his vM) by maximizing an average fi- 
delity. Fiurasek considers a Gaussian distribution of in- 
put coherent states. His average fidelity includes, first. 
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an average over the success probability, normalized to 
an average success probability, averaged over the input 
Gaussian, and, second, an average over the input Gaus- 
sian. He does not quote the probability of success, and 
he only finds the fc = case. He formulates the optimiza- 
tion problem as a semidefinite program, whereas we use 
a simple Lagrange-multiplier maximization. 

It is useful to pause here to summarize properties of the 



fidelity and the success probability. The fidelity (6.231 is 
the probability of the first N + 1 number states in the 
Poisson distribution associated with the coherent state 
\ga). As we anticipated, this fidelity is close to 1 for 
ga <C VN, goes to zero for ga ^ V^, and transitions 
between these two extremes around ga ~ y/N. Indeed, 
we can use the Chernoff bound for the probability in the 
tails of a Poisson distribution with mean g^a^ to bound 
the fidelity in the two extremes [38^ , 



^2-2 



a^ <N : F{a) = l-Pr[r7.>7V + l| g^a^ ] 

2-2\ N+1 



N + 1 

2;r2 ■ 



g^a^ > N : F{a) = Pr[ n < iV | g^a^ ] 



< 



' eg^a^e-3"'^"/^ Y 



(6.24) 



(6.25) 



In this regime all values of k have success probabilities 
that are exponentially small in N, but fc = is the best 
of a sad lot, indicating that it is the best value of k. 
Before investigating the different values of k in detail, 
however, we extend the Kraus operator PjsfKk so that it 
can map outside Sn in a way that increases the fidelity 
and success probability. 

The extension we seek should preserve the phase in- 
sensitivity of Pj^K^ and should not interfere with the 
output of Pj^Kk in the subspace Sn- A glance at 



Eq. (|6.16|) shows that the extension must have the form 

Vn \n) {n + k\. 



the condition 






=Af+l 



/ > T[Tfe 



KlPt^Kk 



E 



Now we impose 



|u„P|n + fc)(n + fc| . 



(6.27) 



-t, 



The term Kj^P^Kk already satisfies the inequality in the 
subspace SN+k spann ed by the first N + k + 1 number 
states [see Eq. (6.19)], and we can maximize the am- 



plifier's success probability by saturating the inequal- 
ity for the second term, i.e., by choosing Vn = 1 for 

n = N + 1,N + 2,... , with the result that 



The width of the transition region can be estimated by re- 
membering that the two-standard-deviation phase-plane 
radius of a coherent state is 1. As a consequence, the 
amplified output begins to contact the number state cut- 
off at N when ga -I- 1 ~ VN and leaves the high-fidelity 
region entirely when ga — 1 ~ yiV. Thus we expect the 
transition from unity fidelity to zero fidelity to occur as 
a varies from (a/ZV — l)/g to {VN + l)/g. 

The fidelity does not depend on k, but the success 
probability does, so the value of k that maximizes the 
success probability can change as a changes. The am- 
plifier map (6.1) cannot depend, of course, on the input 
amplitude, so we must settle on a value of k and ap- 
ply the resulting map to all input coherent states. We 
are most interested in the high-fidelity regime, where the 
leading-order behavior of the success probability (6.22) 



Ti 



p{/\a) = 



iV! 



r,2fc 



[N + ky. g 



2N 



a < y/N/g . (6.26) 



PnKu 



E 

n=N+l 



n) {n + k\ . 



(6.28) 



With this choice, notice that for /c = 0, the additional 
term in T^ is simply the unit operator in the orthocom- 
plement of Sn- 

The extension of the Kraus operator has essentially 
no impact on the operation of the amplifier in the 
high-fidelity input region. It does increase the fidelity 
marginally in the transition region by including in the 
output number-state components with n > N. The 
biggest effect is to increase dramatically the success prob- 
ability in the low-fidelity regime beyond a ~ \/N + k, 
but this improvement is a pyrrhic victory: all it does is 
to allow the amplifier to report that it worked on inputs 
where the output has essentially the same fidelity with 
the target as the input does. 

Using the extended Kraus operators to calculate the 
success probability and the fidelity of the output with 
the target \ga) gives 



a^ -2k I 



Pk{/\a) = {a\TlTk\a)=e-" a 



m 



Fkia) 



(iV + fc)!g2W 



eN{g''a^)+ J2 



{ga\T 



k\a) 



n=N+l 



e-ff " I N\ 1 

Pk{/\a) " pfc(/|a)/e-"^a2fe Iv (iV + A:)!^^ 



a 



2n 



(n + ky. 



cn 



{9'o?)+ E 



n=N+ 



gn^2n 



^ y/n\{n + ky^ 



(6.29) 
(6.30) 
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With the extended Kraus operators, both the success 
probabihty and the fidehty depend on k. In the high- 
fideUty regime, a <C vN/g, the extension terms have 
httle impact: the fidehty hmits to unity, and the suc- 
cess probabihty has the form given in Eq. (6.26), which 
decreases exponentially with N. For a ^ \/N /g, the 
fidelity goes to zero much as it did before. The success 
probability, however, has a new transition that occurs 
at a^ ~ A^ + fc: for cP' ^ N + k, only the extension 
term matters, so the success probability becomes nearly 
the entire probability under a Poisson distribution with 
mean a^, i.e., pfc(/|ck) = Pr[n > iV + fc + 1 | a^ ], and 
this limits to unity as a^ -^ cx). 




FIG. 5: (Color online) Fidelity Fk{a) of Eq. ([63o| (descend- 
ing curves) and success probability pfc(/| a) of Eq. (6.291 (as- 
cending curves) plotted as functions of input amplitude a 
for different extended Kraus operators Tfe with fc = (solid 
lines), 1 (dashed lines), and 2 (dotted lines): (a) g = \/2, 
A/- = 4; (b) p = 3, iV = 9. The inset in (b) illustrates the 
small differences in fidelity, undetectable in the main plot, 
among the three values of k. 



To gain insight into the success probability ( |6.29 l and 
output fidelity (6.30), we plot them in Fig. [S] as a func- 
tion of the input amplitude a for fc = 0, 1, and 2. In 
Fig. ISTa), we take an amplitude gain g — \/2 and iV = 4, 
both of which are too small to see some of the charac- 
teristic features we have discussed. The three fidelity 
curves are approximately unity until a ^ ^/N/g. After 
this point the fidelity decreases to zero. Conversely, the 
three success-probability curves start close to zero and 
rise to unity after \a\ ~ VN/g. Figure p[b) plots the 
same curves for 5 = 3 and A'^ = 9, values big enough to 
see the characteristic features of the two quantities. In 
particular, it is apparent that the fidelity transitions from 
unity fidelity to zero fidelity around a ~ y/N /g = 1, with 



the transition occurring between {VN — l)/g ~ 2/3 and 
{y/N+ l)/g = 4/3, as anticipated. For all three values of 
fc, the success probability in part (c) rises from its initial 
small value to unity, with the rise occurring around the 
second transition at a ~ \/N. 

It turns out that the success probability and fidelity 
for any value of fc are bounded in the following way: 



0<Pfc(/|a)<Po(/|«), 
Fk{a)<Fo{a). 



(6.31) 
(6.32) 



These bounds are illustrated by the examples plotted in 
Fig. [5J and we have proven them analytically. The proof, 
which is tedious, is contained in Appendix[C| The bounds 
confirm that the best value of fc is fc = 0. We also show 
in Appendix [C] that 



Foia) > {ga\PN\ga) ^ 



(6.33) 



which indicates that the fc = extension increases the 
fidelity over that of the restricted Kraus operators. 

If the amplifier quantum operation has Kraus opera- 
tors other than Tq, our analysis shows that these other 
Kraus operators necessarily reduce the fidelity and the 
success probability. This justifies our earlier assumption 
of a single Kraus operator. The best Kraus operator 
is To, and this gives an amplifier quantum operation 
^^ = To tJ,. 

The three plots in Fig. |6j all for fc = 0, have differ- 
ent values of g and N, but roughly the same high-fidelity 
input region: the ratio vN/g = 1 in parts (a) and (c), 
whereas in (b) it is -\/2- The plots include the fidelity 
and success probability coming from the extended Kraus 
operator To and, for comparison, the fidelity and suc- 
cess probability coming from the restricted Kraus oper- 
ator PnKq. Parts (a) and (b) are interesting because 
they have gains typical of that achieved in experiments, 
but the transitions are not very sharp, g and N being 
too small to see the characteristic features of the plotted 
quantities. In part (c), where g = 3 and A^ = 9, the char- 
acteristic features emerge: the extended Kraus operator 
provides a small increase in fidelity through the transi- 
tion region; the success probability using To ascends to 
1 beyond a ~ VN, instead of falling back to nearly zero 
as happens with the success probability that comes from 
PnKq. These plots illustrate the superior qualities of 
the extended Kraus operator Tq; we do not consider the 
restricted Kraus operators again. 

Figure |6] plots two other quantities: the probability- 
fidelity product, po{/\a)FQ{a) — |(ga|To|Q;)p, for our 
phase-insensitive immaculate amplifier and the overlap 



l(«Mh 



-(g-l)^a^ 



The latter can be regarded as the 



fidelity against the target state of a device that does noth- 
ing, i.e., outputs the input. Since nothing can be done 
with unit probability, |(a|(7Q;)p is also the probability- 
fidelity product for a device that does nothing. A mini- 
mal requirement for a useful amplifier is that it be better 
than doing nothing. The plots suggest that, as far as 
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FIG. 6: (Color online) Fidelity Fq(ol) using the extended 
Kraus operator To [Eq. 1 6.30 1] (solid line with fill ed cir - 
cles); corresponding success probability po(/|ci) [Eq. (|6.29l] 



(dashed line with filled squares); fidelity F(a) using the re- 
stricted Kraus operator PnKq [Eq. ( |6.23[ )] (solid descend- 
ing line); corresponding restricted success probability p{/\a) 
[Eq. 6.221 with fc = 0] (dashed line); probability-fidelity prod- 
uct po{7\a)Fo{d) (solid humped line); and overlap |{Q|(;a)|^ 
(dotted line), all plotted as functions of input amplitude a: 
(a) <? = ^/2, iV = 2; (b) 3 = ^/2, iV = 4; (c) <? = 3 TV = 9. 



the probability-fidelity product is concerned, the phase- 
insensitive immaculate amplifier is never better than do- 
ing nothing — indeed, |(ga|To|a)|^ < |(ga|a)|^ foUovirs 
immediately from the fact that Tg is diagonal in the 
number basis with positive eigenvalues bounded above 
by 1 — and approaches that standard only for a > viV, 
where as we have already seen, Tg becomes the identity 
map. For comparison, the probability-fideli ty pro duct for 
an ideal linear amplifier is 1/g^ [see Eq. (2.23)], which 
beats the do-nothing standard for a^ > \ng''/{g — 1)^. 

The key features of the output state of the immacu- 
late amplifier An — Tq Tg are illustrated by the Q- 
distribution plots in Fig. [t] In Fig. JTla), an input state 
within the high-fidelity input region is transformed to an 
output state that is very close to the target output coher- 
ent state. In part (b), however, the input state is beyond 
the high-fidelity input region; the output state gets plas- 
tered against the output arc of radius ViV, producing a 
flattening and distortion along this arc. This distortion 
is very much like that seen in experiments that imple- 
ment immaculate linear amplification "S?! [28l [3Tti33] . (It 
is worth noting that for the uncxtendcd Kraus operator 
PnKq, as a increases beyond ^fN /g, the output state be- 
comes essentially the Fock state |iV).) Parts (c) and (d) 
illustrate the passage through the second transition at 
a ~ \/N ^ as the action of Tg transitions to being that of 
the unit operator. 
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FIG. 7: (Color online) Q distribution of the output state of 
the immaculate linear amplifier given by the extended Kraus 
operator To, with (? = 3 and A'^ = 9, for four amplitudes of 
input coherent state: (a) a — 0.5; (b) a = 1.5; (c) a = 3; 
(d) Q = 5. The (red) dot denotes the center of the input 
coherent state. The transition at input radius y/N/g ~ 1 
is marked by a (red) arc, and its image at the output by 
the (black) arc at radius \/iV = 3. Thus (a) lies within the 
high-fidelity region, and the output looks like an amplified 
coherent state; (b) lies beyond the transition, and its output 
is flattened along the arc of radius \/iV. A second transition 
occurs near a ~ \/iV, as To transitions to being the identity 
operator. Thus (c), lying right in the middle of this second 
transition, has output that is little amplified and is flattened 
along the radial direction, whereas (d), lying well beyond the 
second transition, has output that is nearly identical to the 
input coherent state. 



In Fig. |8]we plot the SNR-based performance measure 
defined in Sec. |IIC[ with the key difference that we have 
two such SNRs: SNRi = (xi) /Sxi = V^a/Sxi for the 
amplitude (radial) quadrature xi and SNR2 = (xi) /Sx2 
for the phase quadrature X2 ((2:2) 



lie 



0). As in Sec. 

the uncertainties in the SNRs are calculated using anti- 
normal ordering, which applies when one intends to mea- 
sure both quadratures [42]. Figurelslplots the SNR quan- 
tities for an input coherent state ja), the target output 
state \ga), and the output of an Tg immaculate ampli- 
fier. As discussed in Sec. II C the right way to take into 
account the success probability of the immaculate am- 
plifier is to multiply the SNRs by the square root of the 
working probability; thus Fig.lSlalso shows plots the root- 
probability-SNRs, v/po(/|a)SNRi and y^po{s^\a)SNR2, 
for the output of the immaculate amplifier. 

Part (a) of Fig. Is] plots these quantities for a gain typi- 
cal of that achieved in experiments. Part (b) has a larger 
gain that shows the characteristic features of these quan- 
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FIG. 8: (Color online) Antinormally ordered quadrature 
SNRs as a function of the input amplitude a. The SNR is de- 
fined as SNRi = (xi) /Sxi for the amplitude (radial) quadra- 
ture xi or as SNR2 = (xi) /Sx2 for the phase quadrature 
X2- Four of the plots are for (i) the input state \a) (dotted 
line) , for w hich SNRi — SNR2 [this is also the bound given in 
Eq. (2.2f 1]; (ii) the output target state \ga) (solid line), for 
which SNRi — SNR2; and (iii) and (iv) the output state of 
the To immaculate amplifier (SNRi: solid line with crosses; 
SNR2: solid line with circles). The other two plots give the 
amplifier SNRs multiplied by the square root of the working 
probability, Y^po(v^|a), as described in the text: the dashed 
line with crosses plots \/po(/JQ)SNRi, and the dashed line 
with circles plots \/po(v^|«)SNR2. For the amplifier plots, 
(a) has g = ^/2, Af = 2, and (b) has ff = 3, iV = 9. 



titles. Within the high-fidehty input region, the output 
SNRs of the amplifier match those of the target output 
state, but they fall away from the target as a moves out 
of the high-fidelity region. The root-probability-SNRs 
show that once the success probability is taken into ac- 
count, the immaculate amplifier does not do as well as the 
input coherent state; it always satisfies the bound (2.21) 



and is not even close to the bound in the high-fidelity 
region. 

One could use other SNR-based performance measures, 
an example being one based on the statistics of number 
of quanta. Doing this can lead to different conclusions. 
In Fig. [9] we consider a number-based SNR defined as 
SNRat — (N) /AN, where A^ = a^a is the number opera- 
tor and AN is the uncertainty in N. Figure |9] shows that 
in terms of SNRjy, first, the output of the immaculate 
amplifier can do better than the target output state and, 
second, the number-based root-probability-SNR, which 
includes the square root of the success probability, can 
exceed that of the input coherent state. The first of 
these improvements seems to arise from the distortion 
of the output state as it leaves the high-fidelity region 
at a ~ vN / g, which is 1 in both plots; this distortion 
amounts to squeezing in the radial direction, as is illus- 
trated in Fig.^b). The second improvement is due to the 
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FIG. 9: (Color online) Number-based SNR measure as a func- 
tion of the input amplitude a. The four plots are SNRjv for 
the input state \a) (dotted line); SNRjv for the output tar- 
get state \ga) (solid line); SNRjv for the output state of the 
To immaculate amplifier (solid line with crosses); and the 
root-probability-SNR measure \/po(v^|q;)SNR]v for the out- 
put state of the Tq immaculate amplifier (dashed line with 
crosses). For the amplifier plots, (a) has g — \/2, N = 2, and 
(b) has 5 = 3, iV = 9. 



same distortion, but is also aided by the increase in suc- 
cess probability, displayed in Fig. ^ for a > y/N. Since 
these improvements arise from effects outside the region 
of high-fidelity immaculate amplification, they seem to 
be incidental to the operation of the device as an immac- 
ulate amplifier. 

We conclude this section by reiterating that in the 
high-fidelity regime, the fc = extended-Kraus-operator 
immaculate amplifier has a success probability [see 
Eq.(|6.26[) with k = 0] 



Po{/\a) 



^2JV 



(6.34) 



This can be regarded as the chief result of this section: 
within the high-fidelity region of operation, an approxi- 
mate phase-insensitive immaculate linear amplifier has a 
success probability that decreases exponentially with the 
size N/g'^ of the high-fidelity input region, with the base 
of the exponential being g^ . This result, for the op- 
timal phase-insensitive immaculate amplifier, indicates 
that the very low success probabilities seen in experi- 
ments j24l [27| |28] , though they might be depressed yet 
further by technical difficulties, arc an unavoidable conse- 
quence of trying to perform phase- insensitive immaculate 
amplification. 
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VII. CONCLUSION 

Imniaculate amplification is an attempt to evade the 
uncertainty principle by praying that the quantum gods 
won't be paying attention all the time. Unfortunately, 
our results indicate that the quantum ones are keenly 
alert and do not suffer hubris gladly. Our chief conclu- 
sion is that immaculate amplifiers, if they operate phase- 
insensitively, cannot achieve both high fidelity to the tar- 
get output state and even reasonably high working prob- 
ability. Indeed, in phase-plane regions where a phase- 
insensitive device amplifies immaculately with high fi- 
delity, the probability that the device works is extremely 
small. The small working probabilities seen in experi- 
ments that implement immaculate amplification are not 
solely a consequence of technical imperfections; they are 
inherent in the nature of phase-insensitive immaculate 
amplification. 

As penance for trying to outfox quantum mechanics, 
we suggest several changes in focus that might reconcile 
the concept of immaculate amplification and quantum 
theory, as well as leading to more positive results than 
those reported here. The first of these is simple: phase- 
insensitive immaculate amplification, with its reduction 
in noise from input to output, might be a step too far; 
perhaps a better sort of device to seek is a probabilistic 
perfect amplifier, which would amplify the symmetrically 
ordered input noise without adding the noise of a (de- 
terministic) ideal linear amplifier. Working probabilities 
for probabilistic perfect amplifiers might be better than 
those we have found for immaculate amplifiers. 

The second change is to abandon hope for invariance 
under phase-plane rotations and of working on more than 
one circle of input coherent states, focusing instead on 
the quite encouraging probabilities we have found for im- 
maculate amplification of sparse collections of coherent 
states on a single input circle. Nondeterministic devices 
have found many uses in quantum information science, 
a notable example being the KLM scheme for linear- 
optical quantum computing (43] . Immaculate amplifiers, 
like the one described formally by Eq. (5.35), which are 



both phase sensitive and amplitude specific, can work 
on sparse collections of coherent states with high success 
probability; they might find application in problems such 
as discrimination of the coherent states used in phase- 
shift keying [TUl E] . There are important questions re- 
garding communications protocols based on such devices: 
How robust are they against amplitude and phase noise in 
the preparation of the input coherent states? How badly 
are rates impacted by the success probability? These 
questions are certainly worth investigating. 

Finally, we suggest a change in the quantum- 
information-science approach to analyzing amplifying 
devices. The literature on immaculate amplification 
has focused on the fidelity of the output with the im- 
maculate target. We have stressed that fidelity can- 
not be considered as a performance measure alone; the 
probability-fidelity product is a better measure of over- 



all performance. Instead of attempting to optimize the 
probability-fidelity product, however, it might be bet- 
ter to develop performance measures suited to specific 
applications. For metrological applications, the root- 
probability-SNR impresses us as an appropriate measure 
of performance. Continuous- variable quantum key distri- 
bution is a communication protocol that might use im- 
maculate amplification and where key rates are an obvi- 
ous performance measure. Some steps have been taken to 
optimize key rates in this context [lH |45] , but more work 
is needed. To paraphrase Emerson, a foolish fidelity to 
fidelity is the hobgoblin of small minds [17] , that is each 
application begs for its own performance measure. 
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Appendix A: Linear dependence of coherent states 
on a circle 



We review the linear dependence of the continuum of 
coherent states on a phase-space circle of radius a cen- 
tered at the origin. The reader should also consult the 
appendix of Ref. [Mj . 

A coherent state is represented in the number basis by 



J<P 



)=e-V2^ 



n=0 



(Al) 



The coherent states \ae^'^), < (j) < 27r, on a circle of 
radius a are complete, but they are not linearly indepen- 
dent. 

These states are linearly dependent, as we can see from 



2ir 



27r 



■e~™'^|ae''^ 



) = 



~aV2_^ 



nl 



0, 



n >0, 
n<0. 



(A2) 

The vanishing of the integral for n < shows the states 
are not linearly independent. 

That these states are complete follows immediately 
from expanding any vector as 



n=0 •' 



J<t>\ 



(A3) 
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where the function x{4') has Fourier representation 



This brings A into the form 



({(p) = Yl ^"^ 



-in(p 



(A4) 



n=0 



with the positive Fourier coefficients uniquely determined 
to be 



(A5) 



Xn = e'^'/2^(n|V'), n>0, 



and the negative Fourier coefficients arbitrary. That the 
negative Fourier coefficients can be changed arbitrarily 
without changing ji/)) expresses the linear dependence of 
the coherent states on a circle. 



Appendix B: Rotationally symmetric quantum 
operations 

The superoperator that effects a rotation by 6 in the 
phase plane is 

n{e) = e'"^^^ 6-''^°*° = Yl e'^""'"^''|n)(n| |m)(m| . 

(Bl) 

A quantum operation A is invariant under rotations if it 
commutes with 7^(6') for ah 9, i.e., TZ{0) oA = Ao TZ{0). 
The symmetry condition implies that A has the form 



A = YJ2^nl\n + k){n\Q\m){r 



(B2) 



k n,m 



That ^ is a quantum operation, i.e., is completely posi- 
tive, implies that A^^' is a positive Hermitian matrix and 
thus can be diagonalized by a unitary matrix: 



A^i=Y.^run^u. 



(k)jj jj{k)* 



(B3) 



yl=^M/'=)0M/'=)^ 

k,l 



(B4) 



where the operators 



M, 



(k) 



EA*'^^n'^i"+^)H 



(B5) 



are orthogonal Kraus operators. Invariance under 
rotations manifests itself as the requirement that these 
Kraus operators have nonzero number-basis matrix 
elements only in one diagonal strip specified by the 
integer k. 



Appendix C: Optimal success probability and fidelity 



In this Appendix we show that the success probabil- 
ities and fidelities of Eqs. (6.29) and (6.30) satisfy the 
bounds (loTl), KT2\, and (fesl. 



We first show the inequalities 

Pk{^\a) >pfe+i(/|a), 
Fk{a) >Fk+i{a); 



(CI) 
(C2) 



from these, we can also conclude that pk{^\oi)Fk{a) > 
Pfc+i(/|a)i^fc+i(a). This proves that the best success 
probability and fidelity are achieved at fc = 0, i.e., by the 
Kraus operator Tg. 



The success-probability inequalities (CI) follow 
straightforwardly from the difference 



■^tx, xt 



fc ^ T[T,. - Tl^,T,+, 



N 



N\ 1 ^{n + k) 



{N + ky.g 



.9^" 1 



n=0 



n 1 



,„2^^ nl - \^ M^,^^.2,\- + k){n + k\>0. (C3) 



N + k+lg^ 



The manifest positivity of Qk means that {a\Qk\a) > 0, which is the inequality (CI). 



To show the fidelity inequalities (C2), we begin by writing Kraus operator (6.28) in the form 



Tffe = X!-^'=('' 



n=0 



{n + ky. 



n){n + k\ 



(C4) 
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where 



iV! 



/fe(n) = < 



nl 



(71 + fc)! 



n = 0, . . . , A^, 

n = A^ + l,7V + 2, 



(C5) 



Notice that fkin) does not decrease with n for n < N , reaches it maximum value ai n — N , and then is a nonincreasing 
function of n for n > N. 
Using fk{n), we can write 



(.galTfcIa) = e-(^'-i)^'/2a'=E[g"/fcH] , 



(C6) 



where E denotes an expectation value with respect to the Poisson distribution Pr[n | a^ ] = |(n|a)p = P„. We also 
have 



p,{/\a) ^ {a\TlT,\a) = a^''E[f,{n)]. 



Thus the fidelity (6.30) can be put in the form 



Fk{a)=e 



-fg--i)|a|- (E[g"/fcH])' 
E[/|(n)] 



(C7) 



(C8) 



For any fc = 0, 1, . . . , we define 



h[n) I n = N + l,N + 2,. 



(C9) 



V"- + ^ + 1 



Notice that hk{n) is a nonincreasing function of n. 
The fidelity inequality (|C2|) equivalent to 



LHS = (E[g"A.+i(n)])'E[/fc2(n)] < (E[5"/fc(n)])'E[A.+i(n)2] = RHS . 
Since, by the Schwarz inequality, 

LHS < E[g"/fe(n)]E[5"/fc(n)/i2(„)]E[/2(n)] ^I, 
we can achieve our objective by showing that T < RHS or, equivalently, that 

ng'Vkin)hl{n)nfi{n)] < E[5"/,(n)]E[/2+i(n)] . 



C30 n—1 



Equation (Cf2| can be written as 

> ^ G{m, fi) = Xm '^("'' ^^ + ^("' "') ' 

?-n,n n— m—0 

where 

G(m,n) = PnPmfk{n)fk{m)hl{m)[g"'fk{n) - g'' fk{m) 



(CIO) 
(Cll) 
(C12) 

(C13) 

(CM) 



In the final form of Eq. ( C13 1, we use the fact that G{n, n) = to exclude the terms along the diagonal from the sum. 



Now what we show is that 

G{m,n) + G{n,m) = P„P„Jfe(n)A.(m)Lg™A(n) - g" fk{m)][hl{m) ~ hl{n)] 



(C15) 



is never positive for n > m. There are three cases to consider. First, when m < n < N, hkim) — hk{n), so the 
quantity (C15) vanishes. Second, when m < N < n, h^.(m) > hk{n) and 



5"/fc(") - 9"fkim) = 5™ ( Mn) ~ p^fk{N) ) < , 



(C16) 
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so the quantity (C15l is not positive. Third, when N < m < n, hk{m) > hk{n) and g"^fk{n) < g^fk{m), so the 



quantity (C15) is not positive. This completes the proof of the inequahties (C2) 



Now we estabhsh the bound (6.33) by writing the fidehty FQ{a) of Eq. (6.30) as 



/ 2-2^ , Af Y^ y " 



I "^^ n -2n \ 



oo _2„ 

/ 2-2\ , 2N V^ O; 

ew(.9 a )+.9 2^ ^T 



n=Af+l 



n=7V+l 



n 



(C17) 



where the first incquahty follows from using g < g g^ in the denominator. This establishes the bound (6.33). 
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